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Abstract. Following Flory's ideality hypothesis the chemical potential of a test chain of length n immersed 
into a dense solution of chemically identical polymers of length distribution P(N) is extensive in n. We 
argue that an additional contribution 5p c (n) ~ +l/py/n arises (p being the monomer density) for all P(N) 
if n <g (N) which can be traced back to the overall incompressibility of the solution leading to a long-range 
repulsion between monomers. Focusing on Flory distributed melts we obtain 8fi c (n) m (1 — 2n/ (N)) j ' p\fn 
for n <JC (N) 2 , hence, 5fi c (n) ~ —1/py/n if n is similar to the typical length of the bath (N). Similar- 
results are obtained for monodisperse solutions. Our perturbation calculations are checked numerically by 
analyzing the annealed length distribution P(N) of linear equilibrium polymers generated by Monte Carlo 
simulation of the bond-fluctuation model. As predicted we find, e.g., the non-exponentiality parameter 
K p = 1 — (N p ) /p\ {N) p to decay as K p » (N) for all moments p of the distribution. 

Key words. Chemical potential - Polymer melts - Equilibrium Polymers 

PACS. 61.25.H- Macromolecular and polymers solutions; polymer melts - 82.35.-x Polymers: properties; 
reactions; polymerization - 05.10.Ln Monte Carlo methods 

1 Introduction cation of this mean-field result is that density fluctuations 

are small beyond the screening length £, hence, negligible 
One of the cornerstones of polymer physics is Flory's ide- g]. The size of a chain segment of arc- length s of a test 
ality hypothesis [11213] which states that polymer chains chain of length n plugged into a melt of chemically identi- 
in the melt follow Gaussian statistics, i.e. they are random cal A^-polymers of (normalized) length distribution P(N) 
walks without long range correlations. The official justifi- 



Wittmer: Non-extensivity of the chemical potential of polymer melts 



and mean length (N) [4] scales, hence, as 

R 2 (s) = (r 2 ) = b 2 s for g < s < n< (N) 2 (1) 

with b denoting the effective bond length and g the number 
of monomers spanning the screening length £ [2J. See Fig.[T] 
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for a sketch of some of the notations used in this paper. Fi g- 1- Sketch of the problem addressed: (a) Two chains repel 
The extensivity of the chemical potential jjt c (n) of the test each other in an incompressible melt since the total density p 



chain with respect to n, 



cannot fluctuate. In d dimensions this constraint leads to an 



entropic "correlation hole potential " u* ~ n/pR d (r 



l-d/2 



Hc(n) = /in for g <C n <C (N) , 



one has to pay for joining both chains with n and R(n) being, 
fj, > being the effective chemical potential per monomer, respectively, the length and the typical size of the test chains 
is yet another well-known consequence of Flory's hypoth- 
esis pp. The upper boundary (iV) 2 indicated in Eq. |lj 
and Eq. ([2]) for later reference is due to the well-known 
swelling of extremely large test chains where the bath acts 
as a good solution [11315] . 

Assuming Eq. (J2J) to hold for the iV-chains of the bath, 
dense grand-canonical "equilibrium polymers" [6|7] are 
thus supposed to be "Flory size distributed" , 



|1I9| . This sets a non-extensive deviation 8p c (n) ~ u* (n) with 
respect to the chemical potential of asymptotically long chains, 
(b) We compute within one-loop perturbation theory the lead- 
ing non-extensive contribution to the mean self-interaction en- 
ergy (u n ) of a test chain with respect to its Gaussian reference 
state, u n being the sum of the effective monomer interactions 
v(r) between pairs of monomers k and I — k + s at distance r. 



P(N) = fie-^ N , 



(3) 



where as elsewhere in this paper, temperature and Boltz- 
mann's constant have been set to unity [5]. Eq. §3§ im- 
plies, of course, that (N p ) = pi/ ' pP for the p-th moment 
of the distribution. Strictly speaking, Eq. applies only 
for g -C N <C (N) 2 , but both limits become irrelevant for 
systems of large mean chain length, fi — > 0, where only 
exponentially few chains are not within these bounds. 

Recently, Flory's hypothesis has been challenged by 
the discovery of long-range intrachain correlations in three 
dimensional melts [5110111112113114115116] and ultrathin 



films |5l9ll7j . The physical mechanism of these correla- 
tions is related to the "correlation hole" of density n/R(n) d 
in d dimensions with R(n) being the typical size of the 
test chain (Fig. [Ha)). Due to the overall incompressibil- 
ity of the melt this is known to set an entropic penalty 
u*{n) w n/pR(n) d which has to be paid if two chains of 
length n are joined together or which is gained if a chain 
is broken into two parts [9]. The same effective repulsion 
acts also between adjacent chain segments of length s and 
this on all scales [11113] . In three dimensions this leads to 
a weak swelling of the chain segments characterized, e.g., 
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by large mean chain length. Eq. (JSj) will allow us to demon- 

R 2 (s) c 

1 —£-L = _JL f or g <^ s <ti n (4) strate Eq. ([6]) numerically from the observed non-exponen- 

ts i/s 

tiality of the length distribution of equilibrium polymer 
melts obtained by means of Monte Carlo simulation of a 
standard lattice model [7119] . 

The one-loop perturbation calculation leading to Eq. © 
is presented in Section [2] where we will also address the 
chemical potential of monodisperse melts. Section [3] out- 
lines the numerical algorithm used for the simulation of 
equilibrium polymers. Our computational results are com- 
pared to theory in Section |4l 



with c s = \J 24 / 'ir 3 / pb 3 denoting the "swelling coefficient" 
[10. 18 13 . Similar corrections with respect to Flory's ide- 
ality hypothesis have been obtained for other intrachain 
conformational properties such as higher moments of the 
segmental size distribution [13j . orientational bond cor- 
relations [10118116] or the single chain structure factor 

mug. 

In this paper we question the validity of Flory's hy- 
pothesis for a central thermodynamic property the chem- 
ical potential fi c (n) of a test chain inserted into a three- 
dimensional melt. Our key claim is that the correlation 
hole potential leads to a deviation 

Sfi c (n) = ji c — [in « u*(n) ~ +l/p^/n for n -C (N) (5) 



2 Perturbation calculation 



2.1 General remarks 



Following Edwards [2] we take as a reference for the per- 
that is non-extensive in chain length and this irrespective turbation calculation a melt of Gaussian chains of effective 
of the distribution P(N) of the bath. Covering a broader bond len g th b - Averages performed over this unperturbed 
n-range we will show explicitly for a melt of quenched reference system are labeled by an index 0. The general 
Flory-size distribution that ^ask ^ s ^° com P u te the ratio Q{n)/Qo(n) of the perturbed 

to the unperturbed single chain partition function 

<$Atc(rc)«-%=(l-2/m) for g < n < (A^) 2 (6) 

where we have defined c M = 3c s /8. This correction implies 
for the annealed length distribution of linear equilibrium 



Qo(n) K h 



n . 

^(n - s) dr G(r,s)v(r) (9) 



polymers that (to leading order) wit h the perturbation potential u n being the sum of the ef- 

„ T r ,„„ fective monomer interactions v(r) of all pairs of monomers 

P(N) « Me -/*JV-«^(A0 (7 ) v > 

of the test chain of length n, and G(r, s) denoting the 
Gaussian propagator for a chain segment of length s [SJ. 
where both the lower (g <C N) and the upper limit (N <C The factor n — s in Eq. ^ counts the number of equivalent 



(8) 



(N) 2 ) of validity become again irrelevant in the limit of monomer pairs separated by an arc- length s. The devia- 
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tion 5p c {n) from Flory's hypothesis is then given by the where we have defined the screening length £ following, 
contribution to (u n ) Q which is non-linear in n. The calcu- e.g., eq 5.38 of Ref. [2]. Interestingly, the number g of 
lation of Eq. @ in d dimensions is most readily performed monomers spanning the blob, i.e., the lower bound of va- 
in Fourier-Laplace space with q being the wavevector con- lidity of various statements made in the Introduction, can 
jugated to the monomer distance r and t the Laplace vari- be determined experimentally or in a computer simula- 
able conjugated to the chain length n. The Laplace trans- tion from the low-wavevector limit of the total monomer 
formed averaged perturbation potential reads structure factor and, due to this operational definition, g is 

/cm sometimes called "dimensionless compressibility" |15I16) . 
dn(u„) e- nt (10) 

j=n Foil) stands for the ideal chain intramolecular structure 

' S "' i l G(q,t)i{q) (11) 



{2ir) d t 2 x ' K ' K ' factor of the given distribution P(N) of the bath. The ef- 

, r i ,2 c j.1 i_- • j fective interaction, Eq. (1121) . depends thus in general on 

where the factor 1/t accounts for the combinatorics and ' ^ ' — r b 

nr ja i i f f \2 i j.\ i 1 1 t~\ ■ T i , the length distribution of the melt the test chain is in- 

G(q,t) = l/((aq) +t) represents the louner-Laplace trans- to 

formed Gaussian propagator G(r, s) 016] with a = b/V2d serted - For Flor y distri butcd melts the structure factor is, 
being a convenient monomeric length. e '^'' Si yen by [12] 

F o(?) = ; i~ (14) 
[aqy + [i 

2.2 Effective interaction potential while for monod i S p erS e melts it reads F (q) = Nf D (y) 

with y = N(aq) 2 and /r>(y) = 2(e _y — 1 + y)/y 2 being 

We have still to specify v(q), the effective interaction be- 

Debye's function [2]. We remind that within the Pade ap- 

tween test chain monomers in reciprocal space. This in- 

proximation for monodisperse chains Eq. (|14[) holds with 

teraction is partially screened by the background of the 

fi replaced by 2/N [2]. 

monomers of the bath. It has been shown by Edwards 

Below we will focus on incompressible solutions (g — 

[2112] that within linear response this corresponds to 

l/vp — » 0) where the effective interaction is given by the 
= ~ + Fo(q)p- (12) inverse structure factor, v(q) w l/pF (q), i.e. we will ig- 

nore local physics on scales smaller than the correlation 

The bare excluded volume v indicated in the first term of 

length £ and assume that both the test chain and the 

Eq. (|12[) characterizes the short-range repulsion between 

ch ains of the bath are larger than g. The effective poten- 
the monomers. Thermodynamic consistency requires |20|12|15] 

tial takes simple forms at low and high wavevectors corre- 

that v is proportional to the inverse of the measured com- 

sponding, respectively, to distances much larger or much 

pressibility of the solution 

smaller than the typical size of the bath chains. In the low- 
v ~ ~ ~ ~2p~ ( a /£) (^^) wavevector limit the potential becomes vq = v(q — > 0) = 
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(N) I (N 2 ^ p for general P(N), i.e. vq — 2p.jp for Flory interested in (possibly diverging) contributions linear in 

distributed and vq = f° r monodisperse melts. Long the length of the test chain or independent of it, we can 

test chains are ruled by vq which acts as a weak repulsive freely subtract linear terms (i.e., terms ~ 1/t 2 in Laplace 

pseudo-potential with associated (bare) Fixman parame- space) or constant terms (i.e., terms ~ 1/t) to regularize 

ter z — vqt/ti. As already recalled in the Introduction, test and to simplify u t . Such a transformation leads to 
chains with n > {(N) 2 / (TV)) 2 w (iV) 2 0] must thus swell f d d q 1 2F _1 - (aq) 2 - t 



and obey excluded volume statistics [l|3j . 

The effective potential of incompressible melts becomes 
scale free for larger wavevectors corresponding to the self- 
similar random walks, 



[2n) d t 2 {aqf+t 2{vp + F - 1 ) 

(16) 

where ". . ." stands for the linear and constant contribu- 
tions we do not compute. Converging now for incompress- 
ible melts for d < 2, the latter reformulation will prove 



(aq) 

for 1 la^r h/NS 1 / 2 (1^) useful below. (See Section [231 for the complete reguliza- 



v{q)^^iovl/q^b{N) k '\ (15) 
2p 



tion of the ultraviolet divergency for incompressible three 
dimensional melts.) 



Applying Eq. (fTo) to incompressible Flory distributed melts 
this leads to 



i.e. the interactions decrease as a power law with distance 
and this irrespective of the length distribution P(N). One 

expects that short test chains with n <C (N) see an inter- 2 4 Flory distributed melts 

action potential of effectively infinite bath chains as de- 
scribed by Eq. (fTS"]) . Please note that Eq. (| 1 5[) lies at the 
heart of the power-law swelling of chain segments, Eq. , 

and related properties alluded to above |10ll2ll3ll6j . Ut = Yp^W J (2^^' ^ + ' ' ' ^ 

= ^Wt 2 -l/t)G(r = 0,t) + ... (18) 
2p 

2.3 Ultraviolet divergency , , , „ ™» „ 

J where we have read Eq. (|17[) as an inverse Courier trans- 

,. „_ TT^~k t form taken at r = 0. Remembering that a factor 1/t in 
Coming back to the computation of Eq. (|11[) one realizes 

, . . , ; . . . , . . , , rr , • t-space stands for an integral f n ds in n-space, the inverse 

that a naive perturbation calculation using the effective JU 

Laplace transform of Ut can be expressed in terms of inte- 
grals of the return probability G(r = 0, s) = (4-7rsa 2 )~ d / 2 . 

We obtain, hence, in n-space 

1 1 / «2-d/2 

1 1 ' l-d/2 " 



interaction given in Eq. (|12p is formally diverging at high 
wavevectors in three dimensions (becoming only regular 
below d = 2) due to the monomer self-interactions which 



should be subtracted. Using Eq. (fT5b instead of Eq. {H} SaJn) = —, — -77^ — -r n ' — «„ 

[a — 2){ATx) a l * pa a \ 2 — d/2 / 

even makes things worse due to an additional divergency (19) 
associated with the self- interactions of the blobs whose where Sp c (n) stands for the non-extensive contribution to 
size was set to zero (g — > 0). However, since we are not (un) - Note that the first term in the brackets scales as 
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n 

Fig. 2. Non-extensive deviation of the chemical potential 
<5/i c («.)/e M as a function of the test chain length n. Flory dis- 
tributed (FD) melts are indicated by solid lines, monodisperse 
(MD) melts by dash-dotted and dashed lines for Pade approx- 
imation and full Debye function, respectively. The bold lines 
refer to test chains in melts of constant (mean) chain length 
with (N) — 100. The deviation changes sign at n ~ (N). The 
thinner lines refer to test chains of same length as the typical 
melt chain, n = (N), where the asymptotic Gaussian behavior 
(Sfi c (n) — 0) is approached systematically from below. 

the correlation hole in d dimension. Its marginal dimen- 
sion is d = 2. The second term characterizes the effective 
two-body interaction of the test chain with itself. As one 
expects pQ, its marginal dimension is d = 4. Although 
Eq. (|19p is formally obtained for d < 2 it applies to higher 
dimensions by analytic continuation. 
In three dimensions Eq. ([IT?)) becomes 

Wn)= (4^^( n " /2 - 2 ^) (20) 
which demonstrates finally the non-extensive correction 
to the ideal polymer chain chemical potential announced 
in Eq. © and sketched for 1/fi = (N) = 100 in Fig. g] 
by the bold solid line. (A slightly different demonstration 



chemical potential of polymer melts 

is given below in Section l2~6l ) As anticipated above, the 
first term in Eq. ([20]) dominates for short test chains. It 
is independent of the polydispersity and scales as the cor- 
relation hole potential, Eq. ([5]). The second term dom- 
inates for large test chains with n > 1/(2/1) becoming 
non-perturbative for n ^>> l//i 2 . Please note that these 
extremely long chains are essentially absent in the Flory 
bath but may be introduced on purpose. Both contribu- 
tions to 5)i c (n) decrease with increasing n. They corre- 
spond to an effective enhancement factor of the partition 
function quite similar to the 5fj, c (n) = — (7 — 1) log(n) in 
the standard excluded volume statistics with 7 w 1.16 > 1 
being the self-avoiding walk susceptibility exponent pQ. 

Interestingly, while 8/j, c (n) decreases at fixed (N) it 
increases as 5fj, c (n) = —c^/y/n for a test chain with n = 
(N), as shown by the thin solid line in Fig.[2j The chemical 
potential of typical chains of the bath approaches thus the 
Gaussian limit from below. 

2.5 Equilibrium polymers 

Flory distributed polymer melts are obtained naturally 
in systems of self-assembled linear equilibrium polymers 
where branching and the formation of closed rings are 
forbidden [6I7I8| . Since the suggested correction S/j, c (N) 
to the ideal chain chemical potential is weak the system 
must remain to leading order Flory distributed and Eq. ([B]) 
should thus hold [21]. Using P(N)/P (N) = Q(N)/Q (N) 
one obtains directly the corrected length distribution for 
equilibrium polymers announced in Eq. |5]). Note that 
Eq. ([5]) is properly normalized, i.e. the prefactor fi of the 
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distribution remains exact if Sfj, c (N) is given by Eq. ([6|). formula, Eq. can be performed leading to 
Since the distribution becomes broader the first moment 1 f d 3 q , 2 



increases slightly at given \x: 

(N) =^ 1 (l + c^). (21) 



Sfi c (n) = — / -j^-^ (exp(~7i(ag) 2 ) - f c (q,n)) 
1 f d 3 q 



rf c (q,n) (25) 



2pJ (2tt) 3 - 

where f c (q,n) = (2 (iV) /F (g) - y) (1 -exp(-xy))/y with 
More generally, one expects for the pth moment x = n/ (N) and y = (N) (aq) 2 . The first term in Eq. ([25]) 



(N*) = ^ (l + ^ (2r(p + 3/2) - r(p + 1/2)) 

(22) 

with r(x) being the Gamma function [22j . The non-expo- 
nentiality parameter K p = 1 — (N p ) /p\ (N) p should thus 
scale as 

K p = WpC^^fjl (23) 

with w p = (r{p + 1/2) + y^ppl - 2r(p + 3/2))/p! being 
a p-dependent geometrical factor. Eq. (|2"3"|) will be tested 
numerically in Section |4j 

2.6 Incompressible melts in three dimensions 



corresponds to the infinite bath chain limit {x —> 0) which 
does not depend on the length distribution P(N). The in- 
tegral over f c (q, n) stands for finite-x corrections for larger 
test chains. 

2.7 Monodisperse melts 

We turn now to incompressible monodisperse melts in 
three dimensions. As already mentioned above, the Debye 
function for monodisperse melts can be approximated by 
the structure factor of Flory distributed melts, Eq. (fT4|) , 
replacing /i by 2/N. It follows thus from Eq. © that 
within Pade approximation we expect 



It is instructive to recover Eq. (|2"01) directly in three dimen- S^ c {n) « Cfl (1 — Ax) (26) 

sions. For that purpose we may subtract 1/ (aq) 2 from the 



propagator G(q, t) in Eq. (TIT))) which amounts to take off 
a linear and a constant contribution. Taking the incom- 
pressible limit this yields 



for n <C ^V 2 with x = n/N. This is indicated by the 
dash-dotted line in Fig. (5) If the test chain and the bath 
chains are of equal length, x = 1, this leads to <5/x c (n) = 
Sc^/y^n, i.e. pL c (n) approaches again its asymptotic limit 



u t = I - 2 I F °^ — ——+ (24) from below (thin dash-dotted line). 

* ' (2tt) 3 t 2 t+(aq) 2 (aq) 2 2p K ' V ' 



for a general structure factor Fo(q). Assuming Eq. (TT4")) for 
Fo(q) we obtain by straightforward integration over mo- 
mentum u t = — (87rpa 3 ) _1 (^ — t)/t 3 ' 2 + . . .. After taking 
the inverse Laplace transform this confirms Eq. (|20p . In- 
terestingly, the inverse Laplace transform of the general 



The calculation of the chemical potential deviations for 
the full Debye function can be performed taking advantage 
of Eq. (|25|) . Specializing the formula to the monodisperse 
case this yields after some simple transformations 



5/i c (n) = ^=(l- I c (x)) (27) 
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where the finite- X correction is expressed by the integral numerically I c (x = 1) = 3.19. We obtain thus 

Ic{x) = V*Jo V )rf S^(n) = -2.19c,/V^ (28) 

Eq. can be evaluated numerically as shown in Fig. [H . . 

as indicated by the thin dashed line. 

where the bold dashed line corresponds to a variation of 
n at constant N = 100 and the thin dashed line to a test 

3 Algorithmic issues 

chain of same length as the chains of the bath, x = 1. 

Both lines are bounded by the predictions for Flory dis- The theoretical predictions derived above should hold in 
tributed melts and the Pade approximation of monodis- any sufficiently dense polymer solution assuming that the 
perse chains. chains are not too short. Since the direct measurement of 
The evaluation of Eq. I[27p deserves some comments, the chemical potential of monodisperse chains (discussed 
The underlined bracket under the integral I c (x) defines a in Section 12. 7p requires a delicate thermodynamic inte- 
slowly varying function of y decreasing from 2 to 1 when gration [23124115] we test the theoretical framework by 
y increases from to infinity. Without this slow factor computing numerically the length distribution P(N) in 
the integral can be scaled: It is proportional to x and systems of annealed equilibrium polymers [3T]. The pre- 
evaluates to I c (x) = 2x in agreement with the correction sented configuration ensembles have been obtained using 
term obtained for Flory distributed polymers, Eq. ([6]). The the well-known "bond fluctuation model" (BFM) [1912517] 
integral is mostly build up by the region xy < A with an — an efficient lattice Monte Carlo scheme where a coarse- 
error ~ 1/ ■s/A. For large x ^> 1 only small y contribute grained monomer occupies 8 lattice sites on a simple cubic 
to the integral; the underlined term in the integral can be lattice (i.e., the volume fraction is 8p) and bonds between 
replaced by 2 and we obtain asymptotically I c {x) = Ax monomers can vary in length and direction. All length 
in agreement with the Pade approximation, Eq. (|26p . If scales are given in units of the lattice constant. Systems 
the first subdominant contribution to the integral is also with an annealed size distribution are obtained by at- 
computed one gets I c (x) — Ax — l.li^fx for large x. For tributing a finite scission energy E to each bond which 
small test chains, the integral provides the first correction has to be paid whenever the bond between two monomers 
I c (x) — 2x, i.e., it vanishes for x — > as already noted, is broken. Standard Metropolis Monte Carlo is used to 
In short, we recover the known asymptotic for short and reversibly break and recombine the chains [7 26J. Branch- 
long test chains but the crossover is very sluggish. The ing and formation of closed rings are forbidden. Only lo- 
simple Pade approximation, Eq. ([26p . is off by f» 40% in cal hopping moves have been used since the breaking and 
the crossover region where x sa 1. Note finally that if the recombination of chains reduce the relaxation times dra- 
test chain is a chain of the bath (x = 1) one evaluates matically compared to monodisperse systems |26] . 
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We only present data for one high density where half 
of the lattice sites are occupied (p = 0.5/8). It has been 
shown 13 15J that for this density we have a dimensionless 
compressibility g — 0.24, i.e. the system may be regarded 
as incompressible on all scales and the lower bound of 
validity of the theory is irrelevant, and a swelling coef- 
ficient c s — 0.41. Hence, c p = 3c s /8 w 0.16 for the only 
parameter of the theory tested here. We use periodic simu- 
lation boxes of linear length L = 256 containing 2 20 ss 10 6 
monomers. The scission energy E has been increased sys- 
tematically up to E = 15 which corresponds to a mean 
chain length (N) w 6100. The configurations used here 
have already been tested and analyzed in previous publica- 
tions discussing the non-ideal behavior of configurational 
intrachain properties |18|ll|12|16j . 

4 Computational results. 

The main panel of Figure[3]presents the normalized length 
distribution P(N) for different scission energies E as in- 
dicated. A nice data collapse is apparently obtained if 
P(N) (N) is plotted as a function of the reduced chain 
length x — N / (N) using the measured mean chain length 
(TV). At first sight, there is no sign of deviation from the 
exponential decay indicated by the solid line. The mean 
chain length itself is given in panel (b) as a function of E 
together with some higher moments (N p ) = ^2 N N p P(N) 
of the distribution. As indicated by the dashed line, we 
find (NP) 1/p ~ exp(£/2) as expected fr om standard lin- 
ear aggregation theory [61718] . The data presented in the 
first two panels of Fig. is thus fully consistent with older 




4 5 6 7 
x=N/<N> 



Fig. 3. Chain length distribution P{N) of linear equilibrium 
polymers for different scission energies E obtained at number 
density p = 0.5/8 of the BFM Monte Carlo algorithm: (a) 
The main panel shows the collapse of the rescaled distribution 
y = P(N) (N) as a function of x = N/ (N) for several E as 
indicated. The exponential decay (solid line) implied by Flory's 
ideality hypothesis is apparently (to first order) consistent with 
our data, (b) First four moments of the distribution vs. E. 
(c) Replot of the data of panel (a) in log-linear coordinates 
focusing on short chains. The data points are systematically 
below the exponential decay (solid line) for i<1. 



computational work [7|26j which has let us to believe that 
Flory's ideality hypothesis holds rigorously. 

Closer inspection of the histograms reveals, however, 
deviations for small x <C 1. As can be seen from panel (c), 
the probability for short chains is reduced with respect 
to the Flory distribution indicated by the solid line. This 
depletion agrees, at least qualitatively, with the predicted 
positive deviation of the chemical potential, Eq. ([5]). 

The curvature of — log(P(7V)), i.e. the non-extensive 
deviation of the chemical potential from Flory's ideality 
hypothesis, is further analyzed in Figure |H Motivated by 
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x=uN 

Fig. 4. Characterization of the deviations of P(N) from 
Eq. 10} using the functional (a) V[P(N)] w 5fi c (N) and 
(b) W[P(N)] w 25fx c (N) - 5ix c (2N) defined in the main text 
which should both vanish for perfectly Flory distributed ideal 
chains as indicated by the horizontal lines in the panels. Data 
points for different scission energies E (as indicated) collapse 
if V[P(iV)]/c MV //I and W[P(iV)]/c M ^(2 - l/y/2) are plotted 
as functions of the reduced chain length x = fiN. For small x 
both functionals decay as \j\fx as shown by the dash-dotted 
lines. The bold lines correspond to the full predictions Eq. (|31[) 
and Eq. J33J for V[P(N)] and W[P(N)], respectively. 

Eq. ([7|), we present in panel (a) the functional 

V[P(N)] = - log(P(A0) -fiN + logO*) (29) 

where the second term takes off the ideal contribution to 
the chemical potential. The last term is due to the nor- 
malization of P(N) and eliminates a trivial vertical shift 
depending on the scission energy E. Consistently with 
Eq. (|21|) , the chemical potential per monomer [i has been 
obtained from the measured mean chain length (N) using 

m = (Ny 1 (l + c^/VWi) ■ (30) 



chemical potential of polymer melts 

Note that fj, and 1 / (N) become numerically indistinguish- 
able for E > 7. If the Gaussian contribution to the chemi- 
cal potential is properly subtracted one expects to obtain 
directly the non-Gaussian deviation to the chemical po- 
tential, Sn c (N) « V[P(N)]. Due to Eq. © the functional 
should thus scale as 

V-[P(A0]/WM«(l-2z)/\/S (31) 

with x = fj,N as indicated by the bold line in the panel. 
This is well born out by the data collapse obtained up to 
x 1=3 5. Obviously, the statistics detoriates for x ^S> 1 for 
all energies due to the exponential cut-off of P(N). Unfor- 
tunately, the statistics of the length histograms decreases 
strongly with E and becomes too low for a meaningful 
comparison for E > 9. It is essentially for this numerical 
reason that we use Eq. (|30[) rather than simple large- E 
limit /i = 1/ (N) since this allows us to add the two his- 
tograms for E = 5 and E = 6 for which high precision data 
is available. Otherwise these energies would deviate from 
Eq. pip for large x due to an insufficient substraction of 
the leading Gaussian contribution to the chemical poten- 
tial. Thus we have used to some extend in panel (a) the 
predicted behavior, Eq. ([S]), presenting strictly speaking a 
(highly non-trivial) self-consistency check of the theory. 

Since the substraction of the large linear Gaussian con- 
tribution is in any case a delicate issue we present in panel 
(b) of Fig. 2]a second functional, 

W[P{N)} = 2V[P(N)} - V[P(2N] 

= log[P(2A0 M /P 2 (A0], (32) 
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where by construction this contribution is eliminated fol- 
lowing a suggestion made recently by Semenov and Johner 
[5] . The normalization factor /i appearing in Eq. (|32p elim- 
inates again a weak vertical scission energy dependence of 
the data. Obviously, W[P(iV)] = for perfectly Flory dis- 
tributed chains. Following Eq. ([7]) one expects W[P(N)] = 
26fj. c (N) - Sfi c (2N) and due to Eq. © 
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with x = fiN. Eq. (J33J) is indicated by the bold line which 
compares again rather well with the presented data. 

The functionals presented in Fig. [3] require histograms 
with very high accuracy. That P(N) is only approximately 
Flory distributed can be more readily seen using the "non- 



Fig. 5. Non-exponentiality parameter K p for different mo- 
ments p as a function of mean chain length (iV). K p is finite 
decreasing with chain length as suggested by Eq. (|23[1 . The ver- 
tical axis has been rescaled successfully using the p-dependent 
weights w p indicated in the inset. 

5 Discussion 



exponentiality parameter" K p = 1 — (N p ) /p\ (N) p which Summary. Challenging Flory's ideality hypothesis, we have 

measures how the moments deviate from the Flory distri- investigated in this study the scaling of the chemical po- 

bution. Obviously, K p = for rigorously Gaussian chains, tential of polymer chains with respect to the length n of 

As stated in Eq. (|2"3")) . we expect the non-exponentiality a tagged test chain plugged into a solution of jV-chains 



parameter to decay as K p = WpC^^fJl w l/p\/ (N), i.e. of a given length distribution P(N) with (N) being the 

as the correlation hole potential of the typical melt chain, typical length of the chains of the bath. By means of 

The main panel of Fig. [5] presents K p /w p c p as a func- one-loop perturbation calculations we have demonstrated 

tion of (N) ps using double- logarithmic axes. The for n <C (N) 2 the existence of a non-extensive deviation 

predicted power-law decay is clearly demonstrated by the Sfi c (n) <C 1 with respect to the Gaussian reference. This 

data. Note that the scaling of the vertical axis with the correction becomes universal for small reduced test chain 

p-dependent geometrical factors w p allows to bring all mo- lengths, x = nj (N) <C 1, scaling as Sfi c (n) ~ ir- 

ments on the same master curve. As can be seen from the respective of the length distribution as suggested by the 

inset of Fig. [5] this scaling is significant since w p varies "correlation hole potential" (Fig. [TJa)). For larger x the 

over nearly a decade between wi = and wg ~ 8.6. correction depends somewhat on P(N), as explicitly dis- 

Deviations from the predicted scaling are visible, not sur- cussed for Flory distributed [Eq. ©] and monodisperse 



prisingly, for small (N) < 10. 



melts [Eq. (|27| ]. but remains generally a monotonously 
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decreasing function of n scaling as methods as proposed in [24] . In particular, this may allow 

r , . 1 t / \ \ *i t / \ k<i „ ,„. N to improve the numerical test of the theory for x ^> 1; due 

dn c (n) w — -={l - I c {x)) with 7 c (x) (34) 



changing sign at x s» 1 (Fig. [2]). For the important limit of 
a test chain of same length as the typical chain of the bath, 
X = 1, the deviation from Flory's hypothesis decreases 
in magnitude with n. For Flory distributed or monodis- 
perse chains I c (x = 1) > 1 and the asymptotic limit, 
^c(n) — > fJ-n, is thus approached from below [Eq. (|28[) ]. 
Note that our predictions are implicit to the theoretical 
framework put forward by Edwards [5] or Schafer [5] , but 
to the best of our knowledge they have not been stated 
explicitly before. 

We have confirmed theory by analyzing in Section 2] 
the length distribution of essentially Flory distributed equi- 
librium polymers obtained for different scission energies 
E by Monte Carlo simulation of the BFM at one melt 
density. Albeit the deviations from Flory's hypothesis are 
small (Fig. [3ja,b)), they can be demonstrated by analyz- 
ing — \og(P(N)) as shown in Fig.[5]or from the scaling of 
the non-exponentiality parameter, K p ~ (N), for all 
moments p sampled (Fig. [5]). We emphasize that the data 
collapse on the theoretical predictions, Eqs. (|31I33I23[) . has 
been achieved without any free adjustable parameter since 



to the exponential cut-off [Eq. ((5J)] this regime has been 
difficult to explore using the equilibrium polymer length 
distribution (Fig. g]) . 

Another interesting testing bed for the proposed corre- 
lation hole effect are polymer melts confined in thin films 
of width H <C R(n) [T7|. A logarithmically decreasing 
non-extensive chemical potential contribution has been 
predicted for these effectively two-dimensional systems [915] . 
The non-exponentiality parameter of equilibrium poly- 
mers confined in thin films should thus decay rather slowly 
with chain length. This is in fact confirmed qualitatively 
by the numerical results presented in Fig. [6] obtained us- 
ing again the BFM algorithm with finite scission energy 
described above. Note that the smallest film width allow- 
ing the overlap of monomers and the crossing of chains 
(H = 4) corresponds to an increase of Ky. by nearly a 
decade for the largest chains we have sampled. The de- 
tailed scaling with H is, however, far from obvious. Larger 
mean chain lengths and better statistics are warranted 
to probe the logarithmic behavior for asymptotically long 
chains predicted by Semenov and Johner [9]. Note that if 
confirmed this prediction should influence the phase dia- 
grams of polymer blends in reduced effective dimensions. 

Outlook. Clearly, the presented study begs for a direct Finally, we would like to point out that the presented 
numerical verification of the suggested non-extensive chem- perturbation calculation for dense polymer chains may 
ical potential for a test chain inserted into a melt of mono- also be of relevance to the chemical potential of dilute 



the coefficient is known. 



disperse chains, Eq. (|26|) . In principle, this should be fea- polymer chains at and around the (9-point which has re- 
sible by thermodynamic integration using multihistogram ceived attention recently [57J . The reason for this connec- 
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< L=512 ► 




<N> 



Fig. 6. Non-exponentiality parameter K 2 = 1- (iV 2 ) /2 (A?) 2 
as a function of the mean chain length (TV) for equilibrium 
polymers confined between two hard walls at distance H as 
indicated. The lateral box size is L — 512, the number den- 
sity p — 0.5/8. In qualitative agreement with theory [9] the 
non-exponentiality parameter is found to decay much more 
slowly with (TV) as in the bulk (squares) , i.e. K2 increases with 
decreasing H. Note that in the large-(iV) limit one expects 
K2 ~ H°/log({N}) as indicated by the dashed line. 

tion is that (taken apart different prefactors) the same 
effective interaction potential, Eq. (fT5|) . enters the per- 
turbation calculation in the low wavevector limit. A non- 
extensive correction Sfi c (n) ~ in three dimensions 
is thus to be expected. 

We thank the Universite de Strasbourg, the CNRS, and the 
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